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Abstract 

In the holographic model of QCD, baryons are chiral solitons sourced by D4 flavor 
instantons in bulk of size with A = g^Nc- Using the ADHM construction we 

explicit the exact two-instanton solution in bulk. We use it to construct the core NN 
potential to order N^/X. The core sources meson fields to order y^Nc/X which are 
shown to contribute to the NN interaction to order Nc/X. In holographic QCD, the 
NN interaction splits into a small core and a large cloud contribution in line with meson 
exchange models. The core part of the interaction is repulsive in the central, spin and 
tensor channels for instantons in the regular gauge. The cloud part of the interaction 
is dominated by omega exchange in the central channel, by pion exchange in the tensor 
channel and by axial-vector exchange in the spin and tensor channels. Vector meson 
exchanges are sub dominant in all channels. 



1 Introduction 



Holographic QCD has provided an insightful look to a number of issues in baryonic physics 
at strong coupling A = g^Nc and large number of colors Nc [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. In 
particular, in [1, 2] baryons are constructed from a five-dimensional Shrodinger-like equation 
whereby the 5th dimension generates mass-like anomalous dimensions through pertinent 
boundary conditions. A number of baryonic properties have followed ranging from baryonic 
spectra to form factors [1, 2]. 

At large Nc baryons are chiral sohtons in QCD. A particularly interesting framework for 
discussing this scenario is the D8-D8 chiral holographic model recently suggested by Sakai 
and Sugimoto [11, 3] (herethrough hQCD). In hQCD D4 static instantons in bulk source 
the chiral sohtons or Skyrmions on the boundary. The instantons have a size of order 1 / -s/A 
and a mass of order Nc\ in units of Mkk, the Kaluza-Klein scale [3]. The static Skyrmion 
is just the instanton holonomy in the 2;-direction, with a larger size of order A° [5]. 

In the past the Skyrmion- Skyrmion interaction was mostly analyzed using the product 
ansatz [12] or some variational techniques [13]. While the product ansatz reveals a pionic 
tail in the spin and tensor channels, it lacks the intermediate range attraction in the scalar 
channel expected from two-pion exchange. In fact the scalar potential to order Nc is found 
to be mostly repulsive, and therefore unsuited for binding nuclear matter at large Nc- The 
core part of the Skyrmion- Skyrmion interaction in the product of two Skyrmions is ansatz 
dependent. In [14] it was shown that the ansatz dependence could be eliminated in the 
two-pion range by adding the pion cloud to the core Skyrmions. In a double expansion using 
large Nc and the pion-range, a scalar attraction was shown to develop in the two-pion range 
in the scalar channel [14]. The expansion gets quickly involved while addressing shorter 
ranges or core interactions. 

In this paper we analyze the two-baryon problem using the D4 two-instanton solution [15] 
to order Nc/X. The ensuing Skyrmion-Skyrmion interaction is essentially that of the two 
cores and the meson cloud composed of (massles) pions and vector mesons. At strong 
coupling, holography fixes the core interactions in a way that the Skyrme model does not. 
Although in QCD very short ranged interactions are controlled by asymptotic freedom, the 
core interactions at intermediate distances maybe still in the realm of strong coupling and 
therefore unamenable to QCD perturbation theory. In this sense, holography will be helpful. 
Also, in holographic QCD the mesonic cloud including pions and vectors is naturally added 
to the core Skyrmions in the framework of semiclassics. These issues will be quantitatively 
addressed in this paper. 
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In section 2, we review the ADHM construction for one and two-instanton following 
on recent work in [15]. In section 3, we show how this construction translates to the one 
and two baryon configuration in holography. In section 4, we construct the bare or core 
Skyrmion-Skyrmion interaction for defensive and combed Skyrmions. We unwind the core 
Skyrmion-Skyrmion interaction at large separations in terms of a dominant Coulomb repul- 
sion in regular gauge. Core issues related to the singular gauge are also discussed. In section 

5, we project the core Skyrmion-Skyrmion contributions onto the core nucleon-nucleon con- 
tributions at large separation using semiclassics in the adiabatic approximation. In section 

6, we include the effects of the mesonic cloud to order N^/X in the Born-Oppenheimer ap- 
proximation. At large separations, the cloud contributions yield a tower of meson exchanges. 
In section 7, we summarize the general structure of the NN potential as a core plus cloud 
contribution in holographic QCD. Our conclusions are in section 8. In Appendix A we detail 
the A; = 1, 2 instantons in the singular gauge. In Appendix B, we revisit the core interaction 
in the singular gauge. In Appendix C, we check our semiclassical cloud calculations in the 
regular gauge, using the strong coupling source theory in the singular gauge. In Appendix 
D we detail our nucleon axial- form factor and the extraction of the axial coupling qa- 



2 YM Instantons from ADHM 

The starting point for baryons in holographic QCD are instantons in flat R\ x Rz- In 
this section we briefly review the ADHM construction [16] for SU(2) Yang-Mills instantons. 
Below SU(2) will be viewed as a flavor group associated to D8-D8 branes. For a thorough 
presentation of the ADHM construction we refer to [17] and references therein. 

In the ADHM construction, all the instanton information is encoded in the matrix-data 
whose elements are quaternions q. The latters are represented as 

q = qMa"" , = (^T^ 1) , (1) 

with M = 1, 2, 3, 4, 1 = l2x2, and the standard Pauli matrices. qM are four real numbers. 
The conjugate (q^) and the modulus {\\q\\) of the quaternion, are defined as 

g^^gM(<7^)^ \\qr^q^q^qq^^\q\t^J2^lt, (2) 

M 

Req = ^qo(T , Imq= = 2^ ftcr , (3) 
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where | g| is the determinant of a matrix q. For clarity, our label conventions are: M, N,P,Q E 
{1, 2, 3, 4}, PL, V, p, (7, e {0, 1, 2, 3}, and i, j, A;, Z e {1, 2, 3} with z = x^. The flavor ,SC/(2) 
group indices are a, 6 e {1, 2, 3}. 

The basic block in the matrix-data is the {l-\-k)xk matrix, A, for the charge k instanton 

A = A + B (g) X , (4) 

where A and B are x-independent {l-\-k) x k quaternionic matrices with information on the 
moduli parameters. We define x — xmc'^ and B(g)x means that each element B is multiplied 
by X. A and B are not arbitrary. They follow from the ADHM constraint 

A^ A = 1 , (5) 

where A^ is the transpose of the quaternionic conjugate of A. f is a k x k invertible 
quaternionic matrix. /"^ (g) 1 means each element is multiplied by 1. The null-space 
of A^ is 2-dimensional since it has 2 fewer rows than columns. The basis vectors for this 
null-space can be assembled into an {1 + k) x 1 quaternionic matrix U 

A^U = , (6) 

where U is normalized as 

U^U = 1 . (7) 
The instanton gauge field A/^ is constructed as 

Am = iU^OmU , (8) 

which yields the field strengths 

Fmn = -2r]aMNU^Hf ® r")Btt/ . (9) 
Self-duahty is exphcit from 't Hooft's self-dual eta symbol 

e„MJV for M, AT = 1, 2, 3 

VaMN = -VaNM = S , r a " 

5aM for AT = 4 
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The action density, trF|^^, can be calculated directly from /, without recourse to the null- 
space U and Fmn [18] 

trF|,^ = tflog|/| , (11) 
where □ = d\j, = 9%d1[, and |/| is the determinant of /. 

2.1 k = 1 instanton 

The k = 1 instanton in the regular gauge is encoded in a quaternionic matrix A 



A = 



A 



At^ (At i-x + X)^) , (12) 



which yields 



n = + {XM - XmY , (13) 



after using (5). Here p (=V^A|^) is the size and {Xm} is the position of the one instanton. 
The field strength is 

Fm« = 'yw^ ((^^^_^'';',^^,), H/t , (14) 
which follows from (9) with 

, P = f"^^Wt, S^fM, (15) 

where p^ = A^A and W e SU{2). The action density follows from (14) or (11) 

95 

trF^^ = log/ = , (16) 

[[xm - XmY + P ) 

which gives the instanton number J d'^xtic = 1 by self duality. The k — 1 instanton 
in the singular gauge is detailed in Appendix A. 
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2.2 k = 2 instanton 

A charge 2 {k = 2) instanton in the regular gauge is encoded in a quaternionic matrix A [15] 

\ 

(17) 



A 



/ Ai A2 

-[x-{X + D)] u 

-[x-{X-D)]J 



u 



where the coordinates xm are defined as x — xmc^ , and the moduli parameters are encoded 
in the free parameters Ai, A2, X,D: \Xi\ = pit are the size parameters, \\\2/ (P1P2) G SU (2) 
is the relative gauge orientation, and X ± D is the location of the constituents, u is not a 
free parameter and will be determined in terms of other moduli parameters by the ADHM 
constraint (5). 

Since we are interested in the relative separation we set X = 0, so that 



A = 



Ai 
D-x 

u 



A2 
u 

-D-x 



At 



'a1 {d-xY 



{-D-xy 



which yields 



A^A 



WXif + \\x - Df + \\uf 
[X\X2 + D^u - uW - {x^u + u^x)] ^ 



AJ A2 + D^u - uW - (a;tit + u^^x) 
\\X2f + \\x + Df + \\uf 



■ (19) 



The ADHM constraint (5) implies that each entry must be proportional to 1. The diagonal 
terms satisfy the constraint. The off-diagonal entries are also proportional to 1 provided 
that u is chosen to be 



u — 



DA 
2Td? 



+ , A = Im(AtAi) = ^(A^Ai - a1A2) , 



(20) 



with 7 an arbitrary real constant. The coordinate u is the inverse of the coordinate D. It 
plays the role of the dual distance. Throughout we follow [15] and choose 7 = for a physical 
identification of the moduli parameters. By that we mean a k — 2 configuration which is the 
closest to the superposition of two instantons in the regular gauge at large separation. In 
Appendix A, we briefiy discuss a minimal k — 2 configuration in the singular gauge. 
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Inserting u into (19) yields 

where we introduced the notation • p for two quaternions q and p 

? • P = ^ QmPm ■ (22) 

M 

Pi — VXi ■ Xi are the size parameters, ±£)m the relative positions of the instantons, and 

2x-u^ [(A2 • D){X, ■ x) - (Ai • D){X2 ■ x) - e^^^«(A2)M(Ai)ArL>pXQ] . (23) 

We made use of the identity 

^P-MN ^ ^PM^N _ ^PN^M _ ^PMNQ^Q ^ 

-MAT ^ 1 (-M^AT _ -AT^M) ^ ^1234 ^ ^ _ (24) 

2.3 Explicit Parametrization 

Without loss of generality, we may choose the moduh parameters to be 

Ai=pi(0,0,0,l) , A2 = P2 (^aSin|^|,cos|^|) , ^=(^^,0,0,0^, (25) 

with a = 1,2,3, I e| = + (^2)^ + (^3)^ and ^„ = The spatial x^ axis is chosen 

as the separation axis of two instantons at large distance d. The flavor orientation angles 
(^a) are relative to the Ai orientation. We assign an SU{2) matrix U to the relative angle 
orientations in flavor space 

U = ^^ e^e-^" e SU{2) , (26) 
which is associated with the orthogonal 50(3) rotation matrix R as 

Rab = ^tr {TaUnU"^) 

= 5„6 COS 2 1 ^ I + 2da9b sin^ | ^ | + eabA sin 2 1 ^ | . (27) 
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For instance Rab reads 



f cos 2^3 sin 2^3 o\ 
-sin 2^3 cos 2^3 
1 



/l ^ 

cos 29i sin 29i 
^0 — sin 2^1 cos 29i J 



(28) 



for 9i = 02 — and 6*2 = 6*3 = respectively. Note the double covering in going from SU(2) 
to S0(3). 

In this coordination for the moduli space, 



A = piIm(A^) = pip2(-^^ar"sin| 9\) , 
DA 



X ■ u — 



iT^K P1P2 . I ^1 ;r la 
9 = — = ^—7- sm 9\ 9aT T°- , 

^ Jnl^l (0,-^3,^2,^1) , 



u — 



P1P2 
Um = — — sm 



d 



(29) 
(30) 

(31) 
(32) 



and the inverse potential / ^ is written as 



/" 



9-+ A B 
B 9+ + A^ 



g± = xl+ {xi± 



d 



( cos I 6*1 + - sin I 6* 



9iz + 92X3 - 93X2 



(33) 

(34) 
(35) 



with p = pi = p2. The action density can be assessed using (11). In terms of this notation, 
for the k = 1 instanton in the regular gauge (13), the logarithmic potential log|/| is 



log/± = -^ogg± , 



(36) 



where the subscript ± refers to the position ^| of the instanton along the Xi axis. For the 
k — 2 instanton (33), we have 



log|/_+| = -\og[{g. + A){g+ + A)-B'] 



(37) 
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2.4 Asymptotics 

To understand in details the structure of the k = 2 instanton it is best to work out its 
asymptotic form for the case d/ 1. For that we use (9) in the special case 



(38) 



with 



/o o\ 

-1 
-1 



(39) 



Below, we will show that the field strength Fj^ sources the pion-nucleon coupling in the 
axial gauge — Q for the quantum fluctuations. The asymptotics is useful for a physical 
identification of the coset parameters. 

Near the singularity center with x — D, (18) approximates to 



At 



'X\ «t 



(40) 



whose null vector U is 



U 



( -±u^ \ 

T^-u (Ay + 2piL)t 



V 



1 



Prom (7) and (20) it follows that 



U 



\ 



1 - (^)^|sin2|^|^„(iTW) 
\^ (f)'sin|^|^„(iTW) J 



+ 



0(f)' 



(41) 



(42) 



We have used the exphcit parametrization (25) and (29). We may expand / near the center 



f\x^D 



9+ 
_A2- 

9-9+ 



XiX2+2x-u I ,c) (ly 
9-9+ ^ \d) 



(43) 
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For X — D, the leading contributions to /n, /12, and /21 are of order l/dP while that of /22 
is of order dP. 

Prom (39) and (42) we have 

a.B|_^(-1.0(9^0(9^).(.l...). (44) 

which yields (38) 

Fi,\x^D = -2U'^B{f®T')B'^U\x^D 

= -2 (/nblr^bi + /i2blr^b2 + /2ibk^bi + /22b5r^b2) 



r 



Thus 



= -2— + 0(d-^) . (45) 

5'+ 



i^:U«D ~ -25'''— (46) 
5'+ 

A rerun of the argument for the center x — —D yields 

Ft,\x^-D ~ -2i?'«^ , (47) 

since W ~ (0,0,1). For asymptotic distances d/p^ 1 the k — 2 configuration splits into 
two independent k — 1 configurations with relative flavor orientation R"^. This separation 
makes explicit the physical interpretation of the coset parameters: p the instanton size, 
d the instanton relative separation, u the inverse or dual separation and R their relative 
orientations asymptotically. 



3 Baryons in hQCD 

Baryons in hQCD are sourccd by instantons in bulk. The induced action by pertinent brane 
embeddings and its instanton content was discussed in [3]. The 5D effective Yang-Mills 
action is the leading terms in the 1/A expansion of the DBI action of the D8 brancs after 
integrating out the 5^. The 5D Chern-Simons action is obtained from the Chern-Simons 
action of the D8 branes by integrating F4 RR flux over the 5"^, which is nothing but Nc- 
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The action reads [11, 3] 



S^Sym + Scs, (48) 



Sym = -^J d^^dz tr ]^K-^l'^3^l, + M^^KS'^ 



(49) 



Scs-^l c.r^\A), (50) 



247r2 



where — 0,1,2,3 are 4D indices and the fifth (internal) coordinate z is dimensionless. 
There are three things which are inherited by the holographic dual gravity theory: Mkk, i^, 
and K. Mkk is the Kaluza-Klein scale and we will set Mkk = 1 as our unit, k and K are 
defined as 

K = AAT^^Y^ = AA^cO , K^l + z'' . (51) 

A is the 5D U{Nf) 1-form gauge field and IT^,^ and are the components of the 2-form 
field strength 3^ — dA — iA/\A. uj^^^^\a) is the Chern-Simons 5- form for the U (Nf) gauge 
field 



u;^^ '"f> (A) = tr (^yir + ^yi'^IJ - —A'j , (52) 

The exact instanton solutions in warped x'^ space are not known. Some generic properties 
of these solutions can be inferred from large A whatever the curvature. Indeed, since k ~ A, 
the instanton solution with unit topological charge that solves the full equations of motion, 
follows from the YM part only in leading order. It has zero size at infinite A. At finite A the 
instanton size is of order 1/\/A. The reason is that while the CS contribution of order A° 
is repulsive and wants the instanton to inflate, the warping in the z-direction of order A° is 
attractive and wants the instanton to deflate in the 2;-direction [2, 3]. 

These observations suggest to use the flat space instanton conflgurations to leading order 
in NcXj with 1/A corrections sought in perturbation theory. The latter is best achieved by 
rescaling the coordinates and the instanton fields as 

Am — X^^'^Am , Ao = Ao , 

3^ MN — MN , 3^0M — X}^^3^ OM ■ (53) 
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The corresponding energy density associated to the action (50) reads [3] 



E = Stt^k 



K 



167r2 



+ — / d^xdz 



-ytr4 + i^tr^-l(aM3or-^ 



(54) 



All quantities are dimensionless in units of Mkk- The U(l) contribution Aq follows from the 
equation of motion [3] 



1 



327r2a 



trF, 



MN 



(55) 



The Ao field can be obtained in closed form using (11), 

1 



Ao 



a log I/I . 



(56) 



According to (53) both the size of the instanton p and the distance d between two instantons 
arc rcscaled, i.e. p = VXp and d = \f\d. While the size p is fixed to po (see below) by the 
energy minimization process, the distance is not. Therefore, when discussing the energy at 
the subleading order, the distance d is always short for ^/\d. It will be recalled whenever 
appropriate. The first term in (54) is Stt^k x instanton number, which is identified with the 
bare soliton mass. The second term (= A£^) is subleading and corresponds to the correction 
to the mass or the interaction energy 



K 

A 

K 



d xdz 



^ 2 ^ 1 ^ 

-^tr ^5 + z'trFl - -{dMAof - 32^^otrF2 



= 5X ^ ^^^^^ ' ^ 



3V 
24 



MN 



□ log I/I log I/I , 



(57) 



where we used the self-duality, trF?- = 2trF^| = |trF^jy, and integrated (9m>1o)^ by P^rt 
so that it can be reduced to the form AotrF|^^. 
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3.1 One baryon 

The one baryon solution is the k = 1 instanton. From (13) it follows that 



r' = p' + xi,, (58) 

for k = 1. Wc have set Xj = by translational symmetry. We have also set X4 = as a 
finite X4 cots energy [3]. Thus 



°^^og/=^^. (60) 



The mass correction AM = AE, reads 



AM(p) = 4 / d^^dl ( ^ + '^^^] J''' (61) 

A V 6 3 207r%2 p2 ^ ■ V ; 

It depends on the size p as plotted in Fig. 1. All integrals in AM are analytical, since 
□ log 1/1 and log |/| are simple. For k — 2 the expressions for AM are more involved and 
require numerical unwinding. As a prelude to these numerics, we have carried the integrals 
in (61) both analytically and numerically as illustrated in Fig.l. 

The one-instanton stabilizes for 



with a mass correction 

AM(po - 9.64) - 0.365 . (64) 

We recall that the physical instanton size po = po/\/X following the unsealing as detailed 
above. 
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AM 




Figure 1: AM/Nc'. solid (exact) and dotted (numerical). 



3.2 Two baryon 

The two baryon solution corresponds to the k — 2 instanton. The corresponding potential 
/ for the k — 2 instanton is given in (33) and yields 



trFi = 



a' log I/I 

= -tflog 



/- , p?plsin^ 1 6'|\ / _ , p?pisin^|6l| 
g.{xM) + ' ' 9+{xm) + - ' ' 



~PiP2 (cos I ^1 + ^sin I 9\ 



(P 
2 
d 



rf2 



6iXo + O2X3 - 63X2 



(65) 



Its leading contribution in (54) is 



2 X Stt^k , 



as expected by self-duality. To order NcX the energy of the 2-baryon system is just 2Mo or 
twice the bare soliton mass. There is complete degeneracy in the moduli parameters d and 
9a- This degeneracy is hfted at order NcX^, which is the next contribution in (54). This will 
be detailed below. 

For two parallel instantons | ^| = and the instanton action density (65) reads 



(66) 



The baryon number distribution in space follows from 



B{x) 



167r2 



(67) 
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Figure 2: Defensive Skyrmions: (a) d = 2, (b) d = v^, (c) d = 1, (d) d = 10 



-5 



which integrates to 2. Since the instanton in bulk is locahzed near z p 1/ \/\, we may 
approximate the integral by the value of the integrand for z ^ 0, or B{x) ~ ii F'^^{z ~ 
0)/167r^. In Fig. 2 we show tvF^^ for z = — and pi = P2 — 9.64 for various separations 
d in the {xi,X2) space for two paralell Skyrmions. The separation is in units of the size 
Po = 9.64. For small separations a narrow Skyrmion develops on top of the broad Skyrmion. 
The configuration is maximally repulsive (defensive Skyrmions). 

A paralell and antiparalell Skyrmion (combed Skyrmions) corresponds to the choice 9i — 
92 — and ^3 = | or | ^| = 7r/2. This is a tt rotation along X3 in the S0(3) notation (27). 
The resulting instanton action density (65) reads 



trF: 



-□Mog 



g-{xM) + 



g+ixu) + 



4^5^ 



(68) 



In Fig. 3 we show the baryon density in the plane (xi, ^2) for various separations in units of 
the instanton size with pi — P2 — 9.64. For large separation two lumps form along the x^ 
axis. For smaller separation the two lumps are seen to form in the orthogonal or X2 direction. 
In between a hollow baryon 2 configuration is seen which is the precursor of the donut seen 
in the baryon number 2 sector of the Skyrme model [19] . The concept of d as a separation 
at small separations is no longer physical given the separation taking place in the transverse 
direction. What is physical is the dual distance u in the transverse plane. 
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(c) (d) 




Figure 3: Combed Skyrmions: (a) d = 2.5, (b) d = 1.7, (c) d = \/2, (d) c? = 1 



For two Skyrmions orthogonal to each other, the choice of angles is = ^2 = 0, = ^. 
The corresponding action density is given by (65) 



trF^, = -tflog 



g-{xM) + 



'pIpI sin^ e. 



d^ 

2x, 



g+{xM) +P2 + ~ 

d^ 



—P1P2 ( cos 63 — sin 63 

d 



(69) 



which is also seen to reduce to (66) and (68) for ^^3 = or tt and 63 = 7i/2 respectively. The 
^3 = ^ is our two orthogonal Skyrmions. This configuration is shown in Fig. (4). For small 
separations a narrow Skyrmion develops on top of a broad one, a situation reminiscent of the 
Defensive Skyrmion configuration above. This situation can be seen in many other relative 
orientations and is somehow generic. 



4 Skyrmion-Skyrmion Interaction 

The Skyrmion-Skyrmion interaction in hQCD is of order Nc/ A and it follows from (57) which 
is the second term in (54). The baryon two minimum energy configuration should follow by 
minimizing this contribution in the 6-dimensional coset space p, d, 9. This will be reported 
elsewhere. Instead, we report on the interaction energy between two Skyrmions versus their 
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Figure 4: Orthogonal Skyrmions: (a) d = 2.5, (b) d = 1.7, (c) d = 1.2, (d) d = 0.6. 

separation for a size fixed in the baryon 1 sector and different relative orientations 9a- In the 
adiabatic quantization scheme, 6a are raised to collective coordinates. They are not fixed 
by minimization. This approach will be subsumed here. We note that the mass shift are of 
order A^^-A*^. 



4.1 General 

Consider the case where 9i — $2 — and ^3 7^ 0, with fixed sizes pi — P2 — Po- Here po 
is the value fixed by minimization in the 1 Skyrmion sector (63). In Fig. (5) we show the 
interaction energy {AE — 2AM) /Nc versus the relative distance d in units of the instanton 
size, where 



6A 



J d^xdz ^2 



24 



□ log I/I log I/I 



(70) 



I/I = g-{xM) + ~ 9+[xm) + ■ 



d^ 



d^ 



pIp^ ( cos 6*3 ^ sin $3 



2xo 



(71) 



The interaction energy is repulsive for all values of 9^. The repulsion decreases for 9s in 
the range — > 7r/2, that is from the defensive to combed configuration. The defensive or 
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AE - 2 AM 




12 3 4 



Figure 5: Skyrmion-Skyrmion interaction in regular gauge. 



AE-2AM AE-2AM 




0.5 1.0 1.5 2.0 2.5 3.0 5 1 1 5 



Figure 6: Skyrmion-Skyrmion interaction: Defensive (left) and Combed (right) 

O3 — IT /2 is still repulsive even for small relative distances, as the two Skyrmions separate 
in the transverse direction. In Fig. (6) we show separatly the interaction energy for the 
defensive configuration (left) and combed configuration (right). The repulsion is seen to 
drop by 3 orders of magnitude. 

The core interaction is modified in the singular gauge as we detail in Appendix A and B. 
In Fig. (7) we show the analogue of Fig. (5) in the singular gauge. The switch from repulsion 
to attraction follows from the switch from repulsive Coulomb (regular gauge) to attractive 
dipole (singular gauge) interactions. The plot is versus d which is the rescaled distance in 
units of the rescaled size p. In the unsealed distance d, the dipole attraction is of order 
Nc/\^ and subleading. 

4.2 Interaction at Large Separation 

To understand the nature of the Skyrmion-Skyrmion interaction to order Nd A as given by 
the classical instantons in bulk, we now detail it for large separations between the instanton 
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Figure 7: Skyrmion-Skyrmion interaction in singular gauge 



cores, i.e. d ^ p but still smaller than the pion range (which is infinite for massless pions). 
We recall that the interaction follows from the subleading contribution in (54), which can 
be split 



c[f] 



D[f] 



N,biC[f]+cD[f]) , 
J d^xdz^DHog\f\ , 

- J d'xdziDHog\f\)^iDHog\f\) , 



with b 



377r2 



6-2167r-^ """"" ~ 24 • 

For large separations between the cores or d ^ p, we have from (37) 



log|/- 



-log 



, , . ,g- + g+ -B"^ 



g-g+ 



logg_ - logg+ - A h 

g-g+ 



g-g+ 



g-g+ 



(72) 
(73) 

(74) 



(75) 



after dropping the A'^ contribution as it is subleading in p/d. We note that after fixing the 
size of the single instanton to po and unsealing the distance d as we indicated above, the 
expansion p/d is an expansion in po/{\/Xd). 

The Skyrmion-Skyrmion core interaction follows from 



V = A^[/_+] - AE[f_] - AE[U] , 



(76) 



after subtraction of the classical self-energies which are of order NcX'^. The C[f] contribution 
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to the interaction reads 



Vc = sin' \e\Vcc, + sin' | 9 \ e}Vcp + sin' | + el)Vcj + cos' | ^1 Vb^ , (77) 

with 

Vco. ^ NJ>t I d?xdz^ □' (3-±^\ , (78) 
(P J \ 9-9+ J 



Vc/3 = N^b^ / d'^xdzz' a' ( J , (79) 

= A^c&^ y c^'5?c^2'i' □' (^i^^ , (80) 

^CT = N^bZ- f d^xdz^ □' ( , (81) 
d' J \9-9+J 

where the cross term in i?' drops by parity and we have rescaled the variable xujd xm- 
Thus g± ^ x' + (xi ± I) ^ +/?/(/'. All integrals are understood in dimensional regularization 
that preserves both gauge and 0(4) symmetry. The results are 

Vca = -Vcp = -Vc-y = N,bCl67r\ Vcs = ■ (82) 

d^ 

The D[f] contribution to the interaction reads 

Vn ~ -2bcN, y ( □' log ^_ ) ^ ( □' log ^+ ) . (83) 

The Coulomb propagator 1/0 = — l/(47r'|x+ — X-|') in 4-dimensions. At large separations 
I x+ — 5^-1 ~ 0? and (83) simplifies to 



VD-12SnHcN,^ 



1 



IGtt^ 



d'^xdz 0^ log g 



2 



277riVc 1 



d2 



(84) 



where the || integrates to the baryon charge 1. Vd captures the Coulomb repulsion between 
two unit baryons in 4 dimensions in the regular gauge. This is not the case in the singular 
as we show in Appendix B. 

We note that after unsealing d — y/Xd, Vd ~ A. In regular gauge, this monopole core 
repulsion is the Coulomb repulsion between 4-dimensional Coulomb charges. We show in 
Appendix C that this is the natural extension of the 3- dimensional omega repulsion at shorter 
distances in holography. The repulsion dominates the many-body problem at finite chemical 
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potential as discussed recently in [4, 5] . Indeed, for baryonic matter at large baryonic density 
ub, the energy is dominated by the Coulomb repulsion (84). The corresponding effective 
interaction is 

Ves = lj dxdy (0+0) {x)Vd{x - y) {ct>+(t>) (y) , (85) 
leading to an energy per volume of order N^vi'^^ /\ as in [5]. 



5 Nucleon-Nucleon Interaction: Core 

At large separation, the nucleon-nucleon core interaction can be readily extracted from the 
Skyrmion-Skyrmion core interaction (77) as it is linear in the 5*0 (3) rotation R. Indeed, 
using the standard decomposition [14] 

R<^ ^-{R^^ + 5''^Rs) , (86) 
3 

with 

Rs = tTR, = Si?"'' - S^hi R , (87) 

the spin Rs and tensor Rt contributions respectively, we may decompose the core potential 
as 

V^Vi + VsRs + V^^R^ . (88) 

The scalar T^, spin Vg and tensor Vt contributions can be unfolded by a pertinent choice of 
orientations of the core Skyrmion-Skyrmion interaction. In general, 

V ^Vi + Vs{Acos^\e\-l) +V^^ (^60''-25"^)sin^|e|+3e"^'=rsin2|^| , (89) 

after using the S0(3) parametrization (27) 

= 5"^cos2|^| +2r^sin2|^| +e"''#sin2|e| . (90) 

The axial symmetry V^(6'i, 6*2, 6*3) = V{9i,9^,92) implies that the tensor components of 
the core satisfy V^"^ = V^^, V^'^ = VS\ and V^^ = V^K Thus, V is reduced to 
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^(^1,^2,^3) = Vi + Vs{4:Cos^\9\-l) + {V^^ -VS''){6el-2)sm''\9\ 

+ (V^^ + V^^)i6§i{d2 + Os)) sin^l ^1) + {V^^ - V^^)(3(^2 + ^^3) sin 2| ^|) 

+ (yS^ + VS^) {662638111^1 6\) + {VS^ - VS^){36ism2\ 6\) . (91) 

In particular, 



y(0, 0, 0) - 14 + 3Vs , V{0, 0, n/2) = Vi-Vs- 2{V^' - V^^) , 
V{7r/2, 0, 0) = V,-Vs + 4{Vr^' - V^^) , 



(92) 

(93) 



so that 



V, = - [V{0, 0, 0) + V{0, 0, n/2) + V{0, n/2, 0) + V{n/2, 0, 0)] , (94) 



V{0, 0,0)-- {V{0, 0, n/2) + V{0, n/2, 0) + V{n/2, 0, 0)) 
1 



- VS' = - [V{n/2, 0, 0) - V{0, 0, n/2)] . 



(95) 
(96) 



Using (78)-(81) we deduce the scalar, spin and tensor core contributions in the form 



Vi = - {3Vca + Vcp + 2Vc^ + Vcs) + Vd = Vd, 

4 

Vs = l (^-Vca - Ivcp - Ivc^ + Vcs^ = , 
V^'-V^^^l{Vcp-Vc^)^0. 



(97) 



The off-diagonal tensor Vr core contribution vanishes. This is clear from (77). Indeed (77) 
can be decomposed as 



V = sin^l^l Vca + sin^|^|^Vc;3 + sin^|^|(^2^ + ^3^)Vc;^ + cos^|^| Vc5 + Vd 
= sin^ I ^1 (Vca + Vc-y) + sin^ | e\ P,{Vcp - Vc^) + cos^ | 6\ Vcs + Vd 
= \ i^Vca + Vcp + 2Vc^ + Vcs) + Vd 

+\ [-Vca - \vc^ - \vc, + Vc^ (4cos^| ^1 - 1) 



22 , 22> 



+-^{Vc0-Vc,){<66l-2) sin^l 



(98) 
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Figure 8: Vi,Vs,Vt in regular gauge 
in agreement with (97). In summary 

Vi^Vd, (99) 

and all others vanish. For general distances, we plot Vi, Vs and Vr with (94)- (96) in Fig. (8) 
in the regular gauge. The relative separation d is in units of the core size p — 9.64. 

In the singular gauge, the Vc core contribution to the nucleon-nucleon interaction remains 
unchanged while the Vd contribution changes. As a result, the spin and tensor channels 
remain the same for both regular and singular gauges. The central or scalar channel Vi — Vd 
changes from repulsive Nc/Xd? (regular) to attractive —Nc/\^d^ (singular) asymptotically. 
The flip is from monopole to dipole as we detail in Appendix B. The short distance repulsion 
in the regular gauge is the 4-dimensional extension of the 3-dimensional omega repulsion. 
In Fig. (9) we show Vi, and Vt with (94)- (96) and (158) in the singular gauge. 

6 Nucleon-Nucleon Interaction: Cloud 

To assess the nucleon-nucleon interaction beyond the core contribution we need to do a 
semiclassical expansion around the k = 2 configuration, thereby including the effects of pions 
and vector mesons as quantum fluctuations around the core. We refer to these contributions 
as the cloud. The semiclassical expansion around the k — 2 configuration parallells entirely 
the same expansion around the k — 1 instanton as detailed in [7]. The expansion will be 
carried out in the axial gauge — for the fluctuations. This gauge has the merit of 
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Figure 9: Vi,Vs,Vt in singular gauge 



exposing explicitly the pion-nucleon coupling. All cloud calculations will be carried with the 
background k — 2 instanton in the regular gauge. Some of the results in the singular gauge 
are reported in Appendix C. 



6.1 Pion 

In the axial gauge — for the fluctuations, the pion coupling to the flavor instanton is 
explicit in bulk. Indeed, following the general expansion in [7] we have for the pion-instanton 
linear coupling 

S^-i^ J d^xdzd, {KF^^C^") , (100) 

with the explicit pion field 

= ^a'^n^o , V'o = - arctan^ , (101) 

and Jt^ — Ak/tt. As noted in [7] all linear meson couplings to the flavor instanton are 
boundary-like owing to the soliton character of the k — 2 instanton. Since KFziipo is odd in 
z, for a static instanton, 

S^K [ d'xF:,Ki^o ^ . (102) 
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Here B — ±Zc refers to boundary of the core when using the non-rigid quantization scheme. 
To avoid double counting, for z < Zc the mesons are excluded in the holographic direction. 
Zc plays the role of the bag radius. It will be reduced to —> at the end of all calculations, 
making the non-rigid quantization constraint point-like. 

The linear pion-2-instanton vertex (102) contributes to the energy through second order 
perturbation. Specifically, 



Ak''K{z^Mzc) 
K^Kiz^fM^^f 



Vu = 0,2 / dxdyF^,{x,z,){d,U{xrd,U{y)'')FSM^c) (103) 



J dxdyFt,{x, z,)did^^^^^FfM ^c) , (104) 



for massless pions. At large separations, the field strength Fj" sphts into two single instantons 
of relative distance d and flavor orientation R. At large relative separation d, (104) simplifies 
to 

Vn«^^JX(0)A,Jr(0), (105) 

with Dij = {3didj — Sij )/d^. The spatial component of the axial vector current J a is unrotated 
while is rotated. Prom Appendix D, its zero momentum limit reads 

J-(0) ^ r^{q = 0) = ~kK{z,)^Po{zc) J dxF^^{x, z,) . (106) 
The projected potential Vu yields 

{sihS2t2 1 Vn I SitiS2t2) = — ^ {Siti\J^{0)\Sih)Dij {S2t2 1 Ja^' (0) 1 52^2) 

1 / \ 2 



167r 



if) ^(3(^l-^)(^2-5)-^l-?2)(Tl-f2) , (107) 



where qa — 32«;7rp^/3 is the axial- vector charge of the nucleon as detailed in Appendix D. 
QA ^ NcX^ in hQCD. 

In the Az — gauge, the linear pion-2-instanton vertex (100) yields a tensor contribution 
to the nucleon-nucleon potential 
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This is in agreement with the pseudo-vector one-pion exchange potential 

^^.n ^ , (109) 

if we identify 

g-KNN QA 



(110) 



This is just the Goldberger-Treiman relation which is also satistified by the holographic 
construction in the — gauge and for massless pions. 

In reaching (107) and the relation (110) there is a subtlety. Indeed in (105) the pion 
propagator Dij is supposed to be longitudinal and the axial vector source transverse, so 
that the contraction vanishes. The subtlety arises from the ambiguity in the axial vector 
source at zero momentum and for massless pions as discussed in Appendix D. The contraction 
is ambiguous through 0/0. The ambiguity is lifted by the order of limits detailed in Appendix 
D, which effectively amounts to a longitudinal component of the axial vector source at zero 
momentum. This result is independently confirmed by using the strong coupling source 
theory discussed in Appendix C. 

Finally, the pion coupling (102) in the axial gauge is pseudoscalar and strong and of 
order y^iVc/A. The reader may object that this conclusion maybe at odd with naive 1/A^c 
power counting whereby the pseudovector coupling is of order \fWf./Nf. ^ with the 

extra l/N^ suppression brought about by the 75 in the nucleon axial vector source [14]. In 
strongly coupled models such as hQCD the nucleon source is of order not 1/Nc, and yet 
chiral symmetry is fully enforced in the nucleon sector. hQCD is a chiral and dynamical 
version of the static Chew model of the A for strong coupling [20]. Also, the reader may 
object that the one-pion iteration which is producing a potential of order N^/X, may cause 
an even stronger correction by double iteration of order (Nc/X)^^"^ and so on. This does not 
happen though, since the direct and crossed diagram to order (Nc/X)^^'^ cancel at strong 
coupling. The same cancellation is at the origin of unitarization in n N ^ n N scattering 
(Bhabha-Heithler mechanism). 



6.2 Axials 

The linear vertex (100) also couples vector and axial vector mesons to the 2-instanton solution 
at the core in bulk. For instance, the axial-vector meson contribution follows from (100) by 
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inserting 



(111) 



so that 



(112) 



The sum over n is subsumed. We have used the fact that KFziip2n is odd in z (axial exchange) 
and that the surface contribution at 2; = oo is zero since 

zu ~ ^{^^ is locahzed in bulk to 

leading order in 1/A. 

In second order perturbation, (112) contributes a static potential 

= 2K''K{z,fi^2n{zc)i^2rn{zc) j dxdyF^^{x,z:)/^^^^^ (113) 

At large separations, the field strength Fj" splits into two single instantons of relative distance 
d and flavor orientation R — RJR2. At large relative separation d, (113) simplifles to ^ 



IGtt 
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167r 



didj \ e 
-A. . J LJ. 
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m 



2n 



-m2nd 



E^x(o)(t 



1 + 



+ 



d 
3 



m2nd ' mi 



-Sij I 1 + 



7^(0) 

dq^d j 



d 



-Jf ^(0) , (114) 



where Ja{0) is defined in (106) and the spatial component of the axial vector current Ja is 
unrotated while is rotated. The projected potential Va yields 



{sitiS2t2\VA\SitiS2t2) 



167r 

n ^ ' 

2 / / \ 2 

9 A \^ li^' 



El \ m2nd 



1 



d ml^d^ 



(cti • (T2) (ti • T2) 



167r 



^ ^0 



+ 



d ' m2nd'^ ml d^ 



_g^^ (ihnV e-"'^-" 

167r ^ V ^0 



(Ji • ci)(a2 • 0?) (fi • fa) 
(Ti • d){a2 ■ d) - ((Ti • (T2) (fi • f2) , (115) 



^For simplicity we often omit \z=zc and = t/Jnizc)- 
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which contributes to the spin Vs,a and tensor part Vt,a of the NN interaction. Specifically, 

VS,A «^ Gl^M^^ ' ^ 5J ^TA,2n^^ , 

n n 

GsA,2n= 7=-— ~ 5'vl"^2n/V« , GrA,2n = QA'm'inNK', (Ho) 

with GsA,2n ~ i/iVc/A and GrA,2n ~ \/iVc/A the spin and tensor couplings of the tower of 
axials to the nucleon. 



6.3 Vectors 

For the vector mesons the time component Fq^ contribution is leading in Nc compared to the 
space component Fiz. This is the opposite of the axial vector contribution. For the SU(2) 
part (rho, rho', ...), we have 

1 r -m2„-i\x-y\ 

z-K ^ J \x -y\ 

1 p — m2n-ld 



n 



where J" = J dx2KKF^Q is the unrotated angular momentum in [7]. We note that 

Z=Zc 

R — RJR2 and that R2^J^ — where I2 is the unrotated isovector charge of the nucleon 
labelled 2. The same holds for label 1. Thus 

^ I Q-m2n-ld 

{sitiS2t2\Vv\sitiS2t2) ~ 2^ 4^2n-l ('^1 ' ^2) , (US) 

n 

which is seen to contribute to the isospin part of the central potential 

p — m2n-ld 

y.-y-^Y.GlvM-^^-r^ (119) 



n 



with G'iy,2n-i = V^2ra-i/2 ~ 1 / \/ N^X- This contribution is subleading in the potential. 
Similarly, the U(l) vector contribution part (omega, omega', ...) reads 

J\f2 . . -m2„-id A7-2 ^ -m2„-id 

- ^E«<-.^^-^E^L-.^. (120) 
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where B = J dx^nKF^o is the baryon number introduced in [7]. This contribution to 

° Z=Zc 

the central potential is leading 



n 



-m.2n-ld 

And ' 



G 



V,2n-l — o ^2n-l , 



(121) 
(122) 



with Gy ~ 

For completeness, we quote the spatial contributions from the vectors, both of which are 
subleading in the potential. The SU{2) vector meson contribution is 

V{, = 2K^K{z,f,p2n-i{z,)ij2ra-i{zc) j dxdyF^,{x,Zc)^Tr^ (123) 

At large separations, the field strength Fj" sphts into two single instantons of relative distance 
d and flavor orientation it! = RjR2- At large relative separation d, (113) simplifies to 



E ^^^(0) (^2-1)^ f-% + (0) 



m 



2n-l 



d 



IGtt 



2 3 , ^ ^ 
1 H ; H ^ ^ I "i" 



m2n-id ml„_id'^ 



-m2n-ld 



'-4^'{Q) , (124) 



where (0) = — {4/3) k,K J dxF^^ and the spatial component of the vector current Jy is 
unrotated while Jy is rotated. The projected potential Vy yields 



{SltlS2t2\Vy\SltlS2t2) 

((Ti • d) {B2 ■ d) (fi • T2) 



Idtt \a m2n-ic( 



IGtt 



E (^2n-l)' 



-m2n-ld 



d 



(cti • d)((T2 • d) - {ax ■ CT2) (fi • f2) , 



(125) 



with Jy (0) = S'y^"'. (125) contributes to both the spin Vg y and tensor part V^ y of the NN 
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interaction, 

^ Q-m2n-ld ^ g-m,2„-irf 

^S,V ~ ^ X/ ^SV,2n-l ^ ' ^T,V ^ ^ X/ <^TF,2n-l ^ , 

n n 

ri — 9v'4^2n-l ^ __ 9vi'2n-\ /■\na\ 

'^SV,2n = , <^TV,2n = ^yj^ " ^ 

The holographic description of the nucleon-nucleon potential is consistent with the meson- 
exchange potentials in nuclear physics. Holography allows a systematic organization of the 
NN potential in the context of semiclassics, with the NN interaction of order Nj A in leading 
order. 



7 Holographic NN potentials 

In general, the NN potential in holography is composed of the core and the cloud contribu- 
tions to order Nc/\. For non-asymptotic distances, both the core and cloud contributions 
have a non-linear dependence on the rotation matrix R{U), making the projection on the 
NN channel involved. Formally, the potential (core plus cloud) can be expanded using the 
irreducible representations of SU(2). Specifically 

oo +j 

V{d, C/) = X E ^jrnid) Di^jU) , (127) 

j=0 m=-j 

where Dl^^,{U) are U- valued Wigner functions. For k = 2 the azimuthal symmetry restrics 
m' = m with Vjm = Vj-m- In particular 

^^-(^) = (2jTI) / ^) ^'rrCmiU) . (128) 

The projection on the NN channel follows by sandwiching (127) between the normalized NN 
states D^/^(l) -D^/^(2). While straightforward, this procedure is involved owing to the 
complicated nature of the k — 2 instanton both in the core and in the cloud on R{U). In 
general 

V{d, U) = Kore(t^, U) + Kloud(c?, U) . (129) 
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VcoTe{d, U) is defined as 



where 



Kore((^, U) = AE[f_+] - AE[f_] - AE[U] , 



K 

6A 



dfxdz ( — 
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□ log I/I log I/I , 



(130) 



log |/_+| = - log [{g_ + A) {g+ + A) - B'] , 
log|/±| = -log^i , 



d 



9± 



p^sin^l^l 



A^ 



a=2,3,4 



B^p^ ( cosl^l + ^sinl^l 



d2 



d 



9iz + 62X3 — 63X2 



ycioud{d, U) is defined as 



V,io.M U) = 2k^K'J2^1 J dxdy 



(131) 
(132) 

(133) 

(134) 



(135) 



f^^{xm; -+) a:^ (f - y) Ft^ivM-, -+) + Fo.(xm; -+) A™(f - ^Fo,(?/m; -+) 



-2i^«(xM; -)A:,^(f - y) f^^1(?/m; -) - 2Fo,{xm; -)A^"(f - y)Fo,(i/M; -) 

where F°'^{xm', — ) and F^^{xm', — H) are the field strengths of the k = 1 and the k = 
instanton respectively, 



2 SU{2) 



(xm; -+) = -25„,C/tB (/_+ ® r") Btf/ , 



(136) 
(137) 
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2,1 i"'\i2 2 2,2,2 

g± = x^+ {xi±-\ +p , x„ = X2 + X3 + X4 , 

, p^sin^l^l ^ 2 / 1^1 2 . 
A = ^ ^2 ' S = pMcos|^| + -sin|^| 



9+ + A, 
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^2 -J 



Sin| e| OaT^T^ 



^ sin| e| ^„tV1 



+ 6'2X3 - 63X2 

[/ = , t/tf/ = 1 . 
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The U(l) fields F^oixM] — h) = dzAo{xM; — h) and Fzo{xm; -) = dzAo{xM; -) follow from 
M^M;-+)^^^alog\f^+\ , lo(xM;-) = ^^nlog|/_| . (138) 
The propagators are defined as 

p-m„\x-y\ -mn\x-y\ 

4:TT\x — y\ 4:Tr\x — y\ 

If we were to saturate (127) by j = 0, 1 which is exact asymptotically as we have shown 
both for the core and cloud, then the projection procedure is much simpler. In particular, 
the NN potential simplifies to 

Vnn = V+ + n-r2 V{- + (tI • a2 (F/ + n-r2 ^5") (140) 
+ (3(<?i ■ d){a2 ■d)-ai- 0^2) {V^ + n-r2 V^) , (141) 

with the core contributions 

nTcore = \ [V{0, 0, 0) + 21/(0, 0, n/2) + V{n/2, 0, 0)] , (142) 



i),core ^ 



V^(0, 0, 0) - ^V{Q, 0, 7r/2) - ^V{n/2, 0, 0) 



(143) 



^TTcore = ' = ^ [V^( ^2, 0, 0) - ^ (0, 0, 7l/2)] , (144) 

as detailed above. The cloud contributions Vi, Vs and Vr remain the same. At large distances 
d the core contribution is dominant and repulsive in the regular gauge (84) 

« , (145) 

and subdominant and attractive in the singular gauge (161) 

' l,core ^4 v "/ 
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The dominant cloud contributions are 



p-m2„-id ]vr jl\F 

Kv " E <^lv,.-.-[^ • G.^,,„., . ^^.„-, ~ sjf . (147) 

n 

Vs,A - GsA,2n^^ , GsA,2n = - ^ ^ , (MS) 

^T,^ - 1. ^TA2n^^ , = ^ ~ y ^ , (149) 

from (108), (116), (121), and (119). To order NJ\ we note that = = 1/+ = 0. 



8 Conclusions 

We have extended the holographic description of the nucleon suggested in [3] to the two 
nucleon problem. In particular, we have shown how the exact k — 2 ADHM instanton 
configuration applies to the NN problem. The NN potential is divided into a short distance 
core contribution and a large distance cloud contribution that is meson mediated. This is 
a first principle description of meson exchange potentials sucessfuUy used for the nucleon- 
nucleon problem in pre-QCD [21]. 

The core contribution in the regular gauge is of order Nc/X. It is Coulomb hke in 
the central channel. Remarkably, the repulsion is 4-dimensional Coulomb, a hallmark of 
holography. This repulsion dominates the high baryon density problem in holography as 
discussed recently in [4, 5]. The dominant Coulomb repulsion is changed to subdominant 
dipole attraction for instantons in the singular gauge. 

We have shown in the context of semiclassics, that the meson-instanton interactions in 
bulk is strong and of order \/Wj\. In the Born-Oppcnheimer approximation they contribute 
to the potentials to order Nc/X. These cloud contributions dominate at large distances. The 
central potential is dominated by a tower of omega exchanges, the tensor potential by a tower 
of pion exchanges, while the spin and tensor potentials are dominated by a tower of axial- 
vector exchanges. The vector exchanges are subdominant at large N^. and strong coupling. 
Holography, fixes the potentials at intermediate and short distances without recourse to 
adhoc form factors [21] or truncation as in the Skyrme model [22]. 

The present work provides a quantitative starting point for an analysis of the NN interac- 
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tion in strong coupling and large Nc- For a realistic comparison with boson exchange models, 
we need to introduce a pion mass. It also offers a systematic framework for discussing the 
deuteron problem, NN form factors and NN-meson and NN-photon emissions in the context 
of holography. We plan to address some of these issues next. 
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A Instantons in Singular gauge 

The k = 1 instanton in the singular gauge follows from (8) through a gauge transformation 
— ^ — ^/l^l which is singular at ^ — x — X — 0. This is achieved through the shift 
U ^ Ug, which amounts in general to the new inverse potential 1// = 1 -\- p^/Clf- The 
corresponding action density is 

t,Fl^ = nMog/ = - l&T^HixM - Xm) . (151) 

[[xm - XmY + p^r 

The instanton in the singular gauge is threaded by an antiinstanton of zero size at its center. 
Its topological charge is strictly speaking zero. It is almost 1 if the center x — X is excluded. 
This point is usually subsumed. Singular instantons have more localized gauge fields than 
regular instantons. 

The ADHM solution for A; = 2 in singular gauge is not known. Following the k — 1 
argument, we may seek it from the regular gauge by applying a doubly singular gauge 
transformation 

g-' = gl^gZ'^i+t , (152) 

which is singular at the centers ^± = a; ± D = in quaternion notations. This amounts to 
shifting U Ug in the ADHM construction. We guess that (152) yields the new inverse 
potential 

^ {\x-D\\x + D\) ' ^^^^^ 
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in the singular gauge. As a result, the instanton topological charge is 



tr Fl^ = log / ^ D^og / - IGtt^ (5(e+M) + (^(^-m)) 



(154) 



The k — 2 ADHM density is now threaded by two singular anti-instantons at ^± — 0. 
Singular instantons have more locahzed gauge fields Am than regular instantons. While 
this point is of relevance for gauge variant quantitities, it is irrelevant for gauge invariant 
quantities with the exception of the topological charge. This point is important for the 
central nucleon-nucleon potential as we now explain. 

B Core in Singular gauge 



In the singular gauge we substitute |/| as (153), which results in 

□ log|/|^nlog|/| + ^ + ^. 



□ ^ log I/I ^ log I/I - IGtT^ (5(e+M) + 5(e-M)) 



(155) 
(156) 



This gives extra contributions in addition to the result in regular gauge after the subtraction 
of the self-energy 



AE 



AE + AE, 

K \ / SV- 
GA 



S 1 

7^2 



24 



d xdz 



327r2{niog|/±| (5e+M) + <^(e-M)) -□log|/+|5(e+M)-niog|/-|5(e-M)} 



+167rM5(e-M)i + 5(e+M)i 
s+ s- 



(157) 



which comes from the second term in (70) while the first term in (70) remains the same. 
/±, /+, and /_ are short for the expressions in the regular gauge in (36)-(37). Thus 



7V,277r 



7V,277r r 4 



d xdz 



+ 



□ log I/I (5(e+M)+5(e- 



16 4 
p2 d> 



d^ 



2d' 



d\2d'' + hd'p' + Ap") + cos' 1^1 (-2rf' - Ad'p' - W cos(2 |^|)) 
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+2dV(d^ + 4rfV + 5p*) sin" 1^1 + rfV(3(^2 + Sp^) sin^ \e\ + 2p^° sin^ |^| 

2' 



d\d^ + p')- cos' 1^1 + d^p'id^ + 2p2) sin' 1^1 + sin^ |^| J | (158) 
iV.27.V + 4cos(2|^|)^^^^_,^ (159) 



2A J d^ 

For large d, the monopole contribution in (158) is cancelled by the monopole contribution 
(84) in the regular gauge. This cancellation leads to a dipole attraction in the singular gauge. 
The net dipole attraction in the singular gauge is best seen by noting that (83) now reads 

Vn^-2bcN,J (tflog(l + p'/ei)) i ( □' log (l + p'/e)) , (160) 

where x± refers to the shifted instanton positions. For large separations d/p':$> 1, the leading 
contribution to Vd is 

~6 ~6 

Vd ~ -7687r' hcN^L^ -SlnN^L , (161) 
d^ d^ 

by repeated use of the 4-dimensional formulae 

° ^ = -4^' S\0 + 2n{2{n + 1) - 4)^ . (162) 

This contribution is of order Nc/X^ following the unsealing of d = y/Xd. (161) is dipole-like 
and attractive. As expected, the threading antiinstanton in the singular gauge cancels the 
leading Nc/X repulsive monopole contribution in 4-dimensional Coulomb's law, resulting in 
the attractive dipole-hke contribution (van der Waals). 



C strong Source Theory 

In this Appendix, we check our cloud calculation in the singular gauge using the strong 
coupling source theory used for small cores in [20, 23] and more recently in holography in [6]. 
This method provides an independent check on our semiclassics in the k = 2 sector. 
The energy in the leading order of A is 

E = Ktr I d'xdz + KFl^ +^ I d'xdz Qi^"^/^^? + KF^^ . (163) 
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where, in the region 1 <^ ^, 



with 



^«-^(G- + G'+), (164) 
« -27rV {{e'^'d+j - S'^d+z) G+ + R"^ (e*^5_,- - S'''d-z) G^) , (165) 
A''^^-27i^p^d+aH+ + R'''d_bH_) , (166) 



oo 

G±^kJ2Mz)M±Z)Y^{\^-X±\) , (167) 

n=l 

oo 



n=0 

b,{z) = ^J-—, (j)n{z) = ^d,Mz) (n= 1,2,3, •••) , (169) 



C.l Pion 

The pion contribution stems from 

= I y d^xdzK (a.A^) (a,yi:) , (171) 

where 

a: «i -27rV + , (172) 

with (f)o{z) only. After subtracting the self-energy the pion interaction energy (Vh) is 



1 5, 



- 2 3% .3 R-[dJ.-f) , (173) 

-* -* J* 

after using 29a = —d±a and dropping the surface terms. X+ — —X^ — | and Zc ~ 0. The 
matrix element of (173) in the 2-nucleon state is 
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(SlS2^1^2|^n|SlS2^1^2) = (3(^1 • d)((T2 • ^) - CTi • CT2) (fi • fs) 



after using (siS2tit2|-R"''|siS2^i^2) = h<^iO'2^i '^2- The last equality follows from the canonical 
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TT A'" pseudovector coupling. Thus 



fg7vNN_Y _ _ fgA 



\ M ) 35 \U 



(175) 



where qa — yKTrp^ obtained in [6], and = 4k /tt. This is just the Goldberger-Treiman 
relation following from the NN interaction using the strongly coupled source approxima- 
tion [6]. As noted in Appendix D, our normalization of the axial- vector current appears to 
be twice the normalization of the same current used in [6] . 

C.2 Omega 

The uj contribution stems from 

^v = lj d'^d^K {9 Jo) (dJo) , (176) 



where 



Ao^-^{G. + G^) . (177) 
After subtracting the self-energy the cu interaction energy (Vy) is 

- (2aA47r)^ \x - X_\\x - ^ ^ 

' / dr{47rr') (179) 



tEA-.^. (180) 
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where we used 



K J dzK{z)d^i^n{z)d^i^rn{z) = ml^_^5nm ■ (181) 

The result is in agreement with (120). At large separations, the strongly coupled source the- 
ory and the semiclassical quantization yields the same results. This outcome is irrespective 
of the use of the singular gauge (strong coupling) or regular gauge (semiclassics) . This a 
consequence of gauge invariance. 

At short distances, gauge invariance is upset by the delta-function singularities present 
in the singular gauge. Indeed, for p <S ^ <S 1 the omega contribution stems from 

Ev^lj d'^dzK (dX) (a,lo) , (182) 



with now K 1 and 



Ao ^ L(G'^^t + G^^^) , (183) 

Zola 

"""" = ■ 

from [6]. After subtracting the self-energy the uj interaction energy iV^) is 



y (2aA47r)2 J ((^^ _ ^/2f + x^)^ {{x, + d/2f + xl)' 



27Nc r z 



2 



dxdz 



27rAcZ2 } ((^^ _ 1/2)2 ^ ((^^ + 1/2)2 ^ 



277riVe 



(185) 



We have rescaled the variable xm — xu/d in the second line and carried numerically the 
integration through 



/ 



^2 TT^ 



dxdz 5 ^ 4.9348 — . (186) 

{{x,-l/2f + xl) {{x, + l/2f + xl) 2 



The omega repulsion at short distance is about y^/2 in (84). The discrepancy may be due to 
the singularities introduced in the singular gauge and/or the approximation in the matching 
region p <^ ^ <^ 1. It is worth noting that the standard omega repulsion at large ditances 
(180) transmutes to a 4-dimensional Coulomb repulsion in holography. 
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D Axial Form Factor 



The effective action for the SU (2)- valued axial vectors to order hP follows from [7] 



3 oo 



6=1 n=l L 

z=B 



z=BJ 



(187) 



The first line is the free action of the massive axial vector meson which gives the meson 
propagator 



Amn,ab , 



x) 



(2^ 



-ipx 



p2 + mln 



(188) 



in Lorentz gauge. The rest are the coupling terms between the source and the instanton: 
the second line is the direct coupling and the last line corresponds to the coupling mediated 
by the SU(2) (a, a', ...) vector meson couplings, 



2n 



;i89) 



which is large and of order 1 / \/h since ■02n ~ V^- When p is set to 1/ \/X after the book- 
keeping noted above, the coupling scales like X^/Nc, or \/iVc in the large Nc limit taken 
first 

The direct coupling drops by the sum rule 



(190) 



cta^ip2n = ipo = - arctau z, 

n=l 

following from closure in curved space 

oo 

5(Z -Z')^Y1 ^^n{z)^n{z')K-''\z') . (191) 
n=l 

in complete analogy with VMD for the pion [11] and the electromagnetic baryon form fac- 
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tor [7]. It follows form (190) that 



n=l 



dz(f)o{z) 



The axial vector contributions at the core sum up to the axial zero mode. 
The iso-axial current is, 



(192) 



(193) 



with 



Ql{y,z) = KK¥',,{y,z) 



(194) 



The static axial-iso- vector form factor follows readily in the form 



J^*(g) = J dxe''^-^ Jf{x) 

= (5- - m I dxe^-' E ^^^V^2n(.) Q?(f , z) 



z=B 



(195) 



and is explicitly transverse for massless pions. The zero momentum limit of the transverse 
momentum projector is ambiguous owing to the divergence of the spatial integrand for g = (). 
We use the rotationally symmetric limit with the convention (^^-^ — q^q') — > 26^^/3. Thus 



4(0) = I / dxf,KF\,{x,z)M^) 



by the sum rule (192). Since the rotated instanton yields 

47r2p' 



z=B 



(196) 



dxFff = 



(197) 



the spatial component of the axial- vector reads 



2^arctan^^^. 



zt 



(198) 
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In the nucleon state 



(s't'lJT^mst) = ^^^^^^^l±|^arctan(..)(a^)-'(rO"' 



9^7 + 



^ IgAi'jT'irY, (199) 
where we used {s't'\R'''\st) = -|((7'')^^'(rO«'. Thus 

= , , , , - arctan(ze) — , (200) 

where the hmit p — > is followed by Zc — > 0. It is 2 times qa — y^Trp^ as quoted in [6]. 
This discrepancy maybe traced back to a factor of 2 discrepancy in the normalization of the 
axial- vector current in [6]. 
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